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Abstract. Controlled generation of entangled states of two quantum bits is a 
fundamental step toward the implementation of a quantum information processor. 
In nano-devices this operation is counteracted by the solid-state environment, 
characterized by broadband and non-monotonic power spectrum often 1/f at low 
frequencies. For single qubit gates, incoherent processes due to fluctuations acting 
on different time scales result in peculiar short- and long-time behaviors. Markovian 
noise originates exponential decay with relaxation and decoherence times, Ti and 
T2, simply related to the symmetry of the qubit- environment coupling Hamiltonian. 
Noise with 1/f power spectrum at low frequencies is instead responsible for defocusing 
processes and algebraic short-times behavior. In this article we identify the relevant 
decoherence times of an entangling operation due to the different decoherence channels 
originated from solid state noise. Entanglement is quantified by the concurrence, which 
we evaluate in analytic form employing a multi-stage approach. "Optimal" operating 
conditions of reduced sensitivity to noise sources are identified. We apply this analysis 
to a superconducting V i — SWAP gate for experimental noise spectra. 
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1. Introduction 

The implementation of a universal two-qubit gate involving an entanglement operation 
on two quantum bits represents a necessary step toward the construction of a scalable 
quantum computer [1]. Intense research on solid state nano-devices during the 
last decade has established the possibility to combine quantum coherent behavior 
with the existing integrated-circuit fabrication technology. In particular, based on 
superconducting technologies, a variety of high-fidelity single qubit gates are nowadays 
available [2l[3ll5, two-qubit logic gates [5l|6] and violations of Bell's inequalities [7] have 
been demonstrated, high-fidelity Bell states generated [8]. The recent demonstrations of 
simple quantum algorithms [9] and three-qubit entanglement [10] are further important 
steps toward a practical quantum computation with superconducting circuits. 

The requirements for building an elementary quantum processor are however quite 
demanding on the efficiency of the protocols. This includes both a severe constraint 
on readout and a sufficient isolation from fluctuations to reduce decoherence effects. 
Solid-state noise sources are often characterized by broad-band and non-monotonic 
power spectrum. Similar noise characteristics have been reported in implementations 
based on Cooper-pair-boxes (CPB) (TT], [12], [T3|, in persistent current [H] and phase 
qubits [TS] [16]. Usually, the spectrum of at least one of the noise sources is 1/f 
at low- frequencies [HI [13 [13 [IH] [20]. At the system's eigen-frequencies instead 
(5 — 15 GHz) indirect measurements indicate white or ohmic spectrum [121 [131 [H]- 
Sometimes spurious resonances of various physical origin have been observed [T5l [T6l [21] . 

At the single-qubit level, the effects of the environmental degrees of freedom 
responsible for the various parts of the spectrum have been clearly identifled leading to 
a convenient classiflcation in terms of quantum noise and adiabatic noise effects [I3l [22] . 
Understanding how these mechanisms affect an entanglement-generating two-qubit gate 
is a relevant issue not yet investigated and it is the subject of the present article. 

The picture for a single qubit can be summarised as follows. Noise at frequencies 
of the order of the system's splittings may induce incoherent energy exchanges between 
qubit and environment {quantum noise). Relaxation processes occur only if the qubit- 
environment interaction induces spin flips in the qubit eigenbasis, i.e. for transverse 
noise. Weakly-coupled Markovian noise can be treated by a Born-Markov master 
equation [23]. It leads to relaxation and decoherence times denoted respectively Ti and 
T2 in Nuclear Magnetic Resonance (NMR) [21]. For transverse noise they are related 
by T2 = 2Ti. Longitudinal noise does not induce spin flips, but it is responsible for 
pure dephasing with a decay-time denoted Tj* [24]. In general, both relaxation and pure 
dephasing processes occur and the resulting decoherence time is T2 = [l/(2Ti)-|-l/T2*]~^. 

Since quantum measurements require averages of measurements runs, the main 
effect of fluctuations with 1/f spectrum is defocusing, similarly to inhomogeneous 
broadening in NMR [21] . Fluctuations with large spectral components at low frequencies 
can be treated as stochastic processes in the adiabatic approximation {adiabatic noise). 
The short-times decay of qubit coherences depends on the symmetry of the qubit- 
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environment coupling Hamiltonian. For transverse noise, the time dependence is 
algebraic oc [1 + at^]~^/^, for longitudinal noise it is exponential quadratic oc exp(— 
("static-path" [22] or "static-noise" [I3] approximation). 

The simultaneous presence of adiabatic and quantum noise can be treated in a 
multi-stage approach [22]. In simplest cases, the effects of the two noise components 
add up independently in the coherences time-dependence. Defocusing is minimized 
when noise is transverse with respect to the qubit Hamiltonian [11]. The qubit is said 
to operate at an "optimal point" characterised by algebraic short-times behavior followed 
by exponential decay on a scale 2Ti. 

In the present article we perform a systematic analysis of the effects and interplay 
of adiabatic and quantum noise on a universal two-qubit gate, extending the multi-stage 
elimination approach introduced in ref. [22]. Understanding these effects is crucial in the 
perspective of implementing solid-state complex architectures. Our system consists of 
two coupled qubits each affected by transverse and longitudinal noise with broad-band 
and non-monotonic spectrum. Such a general situation has not being studied in the 
literature. Previous studies concentrated on harmonic baths with monotonic spectrum 
relying on master equation and/or perturbative Redfield approach [25], or on numerical 
methods [26], or on formal solutions for selected system observables [27] . 

We quantify entanglement via the concurrence [28]. To compare with bit-wise 
measurements, single qubit switching probabilities are also evaluated. Our analysis 
is based on approximate analytic results and exact numerical simulations. Our main 
results are: (i) The identification of characteristic time scales of entanglement decay 
due to adiabatic noise, quantum noise and their interplay; (ii) The characterization 
of relaxation and dephasing for an entanglement operation via the time scales Tr, 
rpswAP^ rpswAP rpswAP*^ polut out the dependence of these scales on the 
symmetry of the Hamiltonian describing the interaction between each qubit and the 
various noise sources; (iii) The demonstration that a universal two-qubit gate can be 
protected against noise by operating at an "optimal coupling", extending the concept 
of single-qubit "optimal point". 

The article is organized as follows. In Section 2 we introduce the Hamiltonian model 
for two-qubit entanglement generation in the presence of independent noise sources 
affecting each unit. In Section 3 the general features of the power spectra of these 
fluctuations, as observed in single qubit experiments, are summarized. The relevant 
dynamical quantities are introduced and the multi-stage approach to eliminate noise 
variables and obtain a reduced description of the two-qubit system is illustrated. In 
Sections 4 and 5 we derive separately the effect of quantum noise within a master 
equation approach and the leading order effect (quasi-static approximation) of adiabatic 
noise. Finally, in Section 6 we discuss their interplay and introduce the relevant time 
scales characterizing loss of coherence and entanglement of a universal two-qubit gate 
in a solid-state environment. Results are summarized in table [5] and in table [61 In 



Appendix A the entanglement generating model is derived for capacitive coupled Cooper 



Pair Boxes (CPBs) including fluctuations of all control parameters. In Appendix B 
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Table 1. Eigenvalues and eigenvectors of Ho expressed in the computational basis 
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the effect of selected impurities strongly coupled to the device is pointed out and 
we speculate on the possibility to extend the "optimal coupling" scheme under these 
conditions. 

2. Universal entangling gate 

Entanglement-generating two qubit gates have been implemented based on different 
coupling strategies. In the standard idea of gate-based quantum computation, the 
coupling between the qubits is switched on for a quantum gate operation and switched 
off after it. The easiest way to realize this scheme is to tune the qubits in resonance 
with each other for efficient coupling and move them out of resonance for decoupling. 
Employing a fixed coupling scheme two-qubit logic gates have been implemented [5] 
and high-fidelity Bell states have been generated in capacitive coupled phase qubits [8] . 
A different idea is to introduce an extra element between the qubits: an adjustable 
coupler, which can turn the coupling on and off [29|, |30]. Alternatively, two-qubit gates 
are generated by applying microwave signals of appropriate frequency, amplitude and 
phase [3Tj . 

The core of the entangling operation of most of the above coupling schemes consists 
of two resonant qubits with a coupling term transverse with respect to the qubits 
quantization axis, as modeled by 

"Ho = - Y (y\z ® I2 - Y II ® 0-2^ + Y (Jxx ® • (1) 

Here Pauli matrices and Iq, is the identity, in qubit-a Hilbert space [a = 1,2). 

In our notation, a^z is the qubit-a quantization axis and we put h = 1. This model 
applies in particular to the fixed, capacitive or inductive, coupling of superconducting 
qubits [32], where individual-qubit control allows an effective switch on/off of the 
interaction. 

Eigenvalues and eigenvectors of eq. ([T]) are reported in tabled! In order to maintain 
the single qubit identities, the coupling strength, Uc, must be one-to-two orders of 
magnitudes smaller than the single qubit level spacing, Q. Thus eigenvalues form a 
doublets structure, as schematically illustrated in figure [H The Hilbert space factorizes 
in two subspaces spanned by {|1), |2)} and {|0), |3)}. A couple of qubits described 
by ([1]) is suitable to demonstrate entanglement generation. The system prepared in the 
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Figure 1. (a) Eigenenergies of the uncoupled resonant qubits; (b) levels in the two- 
qubit Hilbcrt space and logic basis of product states, (c) When the coupling is turned 
on the states | H — ) and | — h) mix and an energy splitting Wc^^ develops between 

the eigenstates {|2),|1)}, spanning the SWAP subspace. Product states | ) and 

I + +) weakly mix and split , with ^3 — wq = 2 i/fi + (wc/S)^. The eigenstates { 1 0) , 1 3) } 
span the Z subspace. Longitudinal noise in the computational basis is responsible for 
inter-doublet relaxation processes (effective transverse inter-doublet) indicated by red 
wavy lines (Eq. (O). Transverse noise in the computational basis originates incoherent 
energy exchanges inside each subspace (effective transverse intra-doublet) indicated by 
blue wavy lines (Eqs. (O, dH)). 



factorized state | H — ) freely evolves to the entangled state \il)e) = [| H — ) — i| — ^)]/V2 in 
a time t^. = it /2uc, realizing a y/i — SWAP operation. The dynamics takes place inside 
the {|1),|2)} subspace, which we name " SWAP-subspace" . The orthogonal subspace 
will be instead named " Z-subspace" . 

Fluctuations of the control parameters used for the manipulation of individual 
qubits couple the circuit to environmental degrees of freedom. We consider the general 
situation where each qubit is affected both by longitudinal noise (coupled to aaz) and 
by transverse noise (coupled to aax), as described by the interaction Hamiltonian 

"Hi = [xi + zi a^] (g) I2 - ^Ii ® [x2 (T2x + ^2 (T2z] ■ (2) 

Here and are collective environmental quantum variables coupled to different 
qubits degrees of freedom. For instance, in the case of two CPB-based qubits at the 
charge optimal point [3], the charge operator is cTq,^., and the Josephson operator is aaz- 
Fluctuations of the gate charge are described by a transverse coupling term, x^aax) 
and noise in the superconducting phase by the longitudinal term, z^ (Jaz [IS] (see 
Appendix A for the derivation) 0. The complete device Hamiltonian reads "Hq+^i+^r, 



where "Hr denotes the free Hamiltonian of all environmental variables. 

In order to identify relaxation and pure dephasing processes for the coupled qubit 
setup we project 'Hq + 'Hi in the 4-dim Hilbert space generated by the eigenstates of Tio, 

I In our notation the qubit's quantization axis is aaz, irrespective of the working point. Thus, Xa and 
Za describe physically different processes only at selected operating points. Usually, this is the case at 
the single qubit's optimal points, see [Appendix A] 
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{\i)}, i = 0, 1, 2, 3, where we may rewrite 
Tio = ^uji\i){i 
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where A± = [cos(v2/2) ± sin((^/2)]/-\/2, with tarup = —uJc/{2^l). In (jlj) we distinguish 
"effective longitudinal" and "effective transverse" terms. The first ones are diagonal 
in the eigenbasis {\i)} and are responsible for pure dephasing processes. "Effective 
transverse" terms instead are off-diagonal and originate both intra- and inter-doublet 
relaxation processes. Specifically we have: 

SWAP subspace {|1), |2)}: Longitudinal noise in the computational basis oc (Jaz^a, 
originates "effective transverse" noise in the SWAP-subspace, i.e. the restriction of 
Ho + Hi to {|1), |2)} reads 



2^ ^ 



Z2) 



|1)(2| + |2)(1| 



(5) 



Note that if both qubits were affected by the same longitudinal noise no effective 
transverse noise in this subspace would be present. This situation may occur in the 
presence of totally correlated noise affecting both qubits [33] . 

Z subspace {|0), |3)}: Longitudinal noise in the computational basis originates both 
"effective transverse" and "effective longitudinal" noise in the Z-subspace, i.e. the 
projection of T-Lq + 'Hi on {|0), |3)} reads 

v/^^^ + W2F[|3)(3|-|0)(0|] 



yproj 



+ Ul + ^2 



cos^(|0)(0|-|3)(3|) + sin(p(|0)(3| + |3)(0|) 



(6) 



Effective longitudinal and transverse components are modulated via the mixing angle, 
(f, similarly to a single qubit with operating point ip. 

Inter-doublet processes: The only effect of transverse noise in the computational 
basis is to mix the two subspaces via "effective transverse" inter-doublet terms (FigJT] 
(c)) 



H 



inter 



2(^1+^2) 



A_|2)(0| + A+|2)(3|+h.c. 



+ 2(^1-^2) -A+|l)(0| + A_|l)(3|+h.c. 



(7) 



These inter-doublet terms are responsible, in particular, for relaxation processes from the 
SWAP-subspace to the ground state. We will demonstrate that the resulting "global" 
relaxation time sets the upper limit to all other gate operation times, including other 
decoherence time scales. 
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3. Multi-scale approach to broad-band noise 

The considered entanglement generating operation takes place in the presence of broad- 
band and non-monotonic noise. In this Section we review the multi-stage elimination 
approach to deal with this problem. The method has been introduced for a single 
qubit in ref. [22] where the various approximations have been checked by comparing 
with the exact numerical solution of the system evolution. This approach allowed to 
accurately explain the observed dynamics in different experiments [T3l [H] , confirming 
its appropriateness to deal with the more complex system studied in the present article. 
The method has been extended to a multi-qubit gate in ref. [M], here we summarize the 
main steps. 

The multi-stage elimination approach is based on a classification of the noise sources 
according to their effects and circumvents the problem of a microscopic description 
of noise sources, which are often non Gaussian and non Markovian [351 IM]- In this 
perspective, the required statistical information on the environment depends on the 
specific quantum operation performed and on the measurement protocol. Even if the 
statistical characterization of the environment requires going beyond the second order 
cumulant, often knowledge of the power spectrum of the bath variables, here Xa and Za 
denoted generically as Ea, 

1 r+oo 

SeM = - die--* [{{E^rnM)) + ((^.(O)^.(t)))] , (8) 

is sufficient. Where ((...)) denotes the equilibrium average with respect to "Hr and we 
assumed stationary processes with {{Ea}) = 0. The typical power spectrum reported 
in various single qubit experiments is sketched in figure [21 To make explicit reference 
to some practical situations, in table [21 we summarize the characteristics of transverse 
and longitudinal noise spectra at low- and high-frequencies reported in a CPB-based 
circuit [13] and in the recent experiment on a flux qubit [H]. 

The low-frequency part of the spectra of each variable is 1//, S^^{u) = Ae^/oj. 
The amplitude Ae^ can be estimated from spectral measurements. If 7^ and •jm 
denote respectively the low- and high-frequency cut-offs of the 1/f region, then 

= 7rc^|Jln(7Af/7m)]"\ where (t|^ is the variance, = J^^' ^ S^^f (u) . It can 
be approximated as (t|;^ = J^l'^ ^S]^^{u!) where tm is the overall acquisition time for 
a single data point which results from averaging over several measurement trials. [§| 

The high-frequency power spectrum is usually inferred indirectly from measure- 
ments of the qubit relaxation times under various protocols [121 ESI HI]- From the 
resulting figures the expected spectra for the corresponding quantum variables Ea at 
the relevant frequencies, Q/2tt and uJc/2'k are derived. Experiments tuning the single 
qubit level spacing Q reveal either ohmic [121 E] or white [13] power spectrum in the 

§ The intrinsic high-frequency cut-off of the 1/f spectrum depends on the specific microscopic source 
and it is usuaUy not detectable in experiments. In the present article we discuss measurements protocols 
where the details of the behavior of the power spectrum below 7™ and around jm sltc not relevant and 
results depend logarithmically on the ratio ^Ai/jm- 
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Figure 2. Sketch of the typical power spectrum of the environmental variable Ea 
(logarithmic scale). Measurements of 1/f noise usually extend between 1 Hz and 
0.1—1 MHz, whereas the ohmic or white spectrum region typically ranges around 
5 — 20 GHz [m [M] . The region classified as adiabatic noise and quantum noise are 
indicated. 



Table 2. Characteristics of transverse and longitudinal noise at low- and high- 
frequencies inferred from data reported in ref. [3] for a charge-phase qubit at its double 
optimal point with « 27r x 16GHz (charge noise is transverse and phase noise is 
longitudinal), and in ref. [13] for a flux qubit at the optimal point with i7 « 27r x 5GHz 
(flux noise is transverse and critical current noise is longitudinal) . For the charge-phase 
qubit we considered jM/jm = 10^, results logarithmically depend on this ratio. 





Charge-Phase Qubit 


Flux Qubit 




CTx « 2 X 10^2 ^ 


CTj; w 2 X IQ-'^n 




cr^ « 10-*^ n 




Sx 


s,:{n) w 4 X lo^s-i 


S^{n) w 2 X lO^s-i 




S,{Vt) « 10^ s-1 





GHz range. Evidence of spurious resonances in the spectrum have often been reported, 
they show up as beatings in time resolved measurements [31 HE] . 

In our analysis noise sources belong to three classes. Low frequency noise with 
1/f spectrum is adiabatic since it does not induce transitions, but mainly defocuses 
the signal, we classify it as adiabatic noise. Noise at frequencies of the order of the 
qubits splittings is responsible for dissipation and ultimately for spontaneous decay, 
thus it is classified as quantum noise. Possible resonances in the spectrum pertain to 
the class named strongly coupled noise. Noise sources belonging to different classes act 
on different frequency scales and are treated via specific approximation schemes. The 
distiction can be illustrated as follows. We are interested to a reduced description of 
the 2-qubit system, expressed by the reduced density matrix (RDM), p(t) obtained by 
tracing out environmental degrees of freedom from the total density matrix p^'^'^'^it), 
which depends on quantum (Q), adiabatic (A) and strongly coupled (SC) bath variables. 
Since bath's degrees of freedom belonging to different classes of noise act of different 
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time scales we separate in each part of the interaction Hamiltonian, aaiE^ {i = x,z), 
the contribution from various noise classes as follows 

aaiEa a^.E'^ + a^iE^ + ac^^El^ . (9) 

Adiabatic noise, E^, is typically correlated on a time scale much longer than the inverse 
of the qubit's frequencies, ^q,, then in the spirit of the Born-Oppenheimer approximation 
it can be seen as a classical stochastic field {Ea{t)} = E{t). This approach is valid 
when the contribution of adiabatic noise to spontaneous decay is negligible, a necessary 
condition being t <^ Tf oc SEj^{^a)~^ ■ This condition is usually satisfied at short enough 
times, since 

^Ea 

(oj) is substantially different from zero only at frequencies u <^ Qa- This 
fact suggests how to trace-out different noise classes in the appropriate order. The total 
density matrix parametrically depends on the specific realization of the slow random 
drives E{t) and may be written as pQ'^' 

sc(i) = pQ'SC(t|^(t)). The first step is to trace 
out quantum noise. In the simplest cases this requires solving a master equation. In 
a second stage, the average over all the realizations of the stochastic processes, E(t), 
is performed. This leads to a reduced density matrix for the 2-qubit system plus the 
strongly coupled degrees of freedom. These latter have to be traced out in a final stage 
by solving the Heisenberg equations of motion, or by approaches suitable to the specific 
microscopic Hamiltonian or interaction. For instance, the dynamics may be solved 
exactly for some special quantum impurity models at pure dephasing, when impurities 
are longitudinally coupled to each qubit [35], [36] . The multi-stage elimination can be 
formally written as 

Pit) = Trsc { J V[Eit)] Pim] Tr^ [p^'^^ (t|E(t)) ] } , (10) 

where Trq and Trsc indicate respectively the trace over the Q and SC degrees of 
freedom. In the following Sections we apply the multi-stage approach to the two-qubit 



gate by tracing out first quantum noise and secondly the adiabatic noise. In Appendix B 
the effect of a SC impurity will be analyzed. 



3.1. Relevant dynamical quantities 

We focus on the Vi - SWAP operation | H — ) -> lipe) = [I H — ) - ^^1 — l~)]/v^ which 
generates by free evolution an entangled state at te = t\ jluJc- As a unambiguous test of 
entanglement generation and its degradation due to noise, we calculate the evolution of 
concurrence during the gate operation. Introduced in ref. [28], the concurrence quantifies 
the entanglement of a pair of qubits, being C = for separable states and C = 1 for 
maximally entangled states. For the situations discussed in the present article and 
specified in the following Sections, the two-qubits RDM takes the "X-form", i.e. the 
RDM expressed in the eigenstates basis is non-vanishing only along the diagonal and 
anti-diagonal at any time. Under these conditions, the concurrence can be evaluated in 
analytic form [37] . In general, it depends both on diagonal and off-diagonal elements of 
the RDM. 
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In order to directly compare with experiments where bit-wise readout is performed, 
we also evaluate the qubit 1 switching probability P^wiit), i-e. the probability that it 
will pass to the state |— )i starting from the state and the probability P2(t) of 

finding the qubit 2 in the initial state |— )2- In terms of the two qubit RDM in the 
eigenstate basis they read (Trjp(t) denotes partial trace over qubit i of the two-qubit 
density matrix) 

Pswi(t) = i(-|Tr2p(t)|-)i = i [pii(t) + P22(t)] + PooW 

+ [PssW - PooW] sin^ I + Re[pi2(t)] + Re[po3(t)] siny^ (11) 

P2{t) = 2(-|Trip(t)|-)2 = i [pii(t) +P22(t)] +PooW 

+ [PasW - PooW] sin^ | - Re[pi2(t)] + Re[po3(t)] sin^. (12) 

For preparation in the state | H — ), in the absence of external fluctuations the above 
probabilities read 

= . P,(t) = i±^. (13) 

The cyclic anti-correlation of the probabilities signals the formation of the entangled 
state, as reported in various recent experiments O [6l [H 138] . 

Both the concurrence and the switching probabilities depend on combinations 
of populations and coherences in the eigenbasis. Therefore the relevant time scale 
to quantify the "quality factor" or the efficiency of the universal two-qubit gate is 
not simply related to a specific RDM element, as a difference with a single qubit 
gate, where the decay time of the qubit coherence quantifies the quality factor of the 
operation (with T2 due partly to relaxation processes (2Ti), partly to Markovian pure 
dephasing processes T2 or originated from inhomogeneous broadening). In the following 
we will study the time dependence of coherences and populations, and identify the 
environmental processes (transverse, longitudinal, low frequency, high frequency, etc.) 
which originate various decay times. Based on this analysis we will discuss the resulting 
effect on the decay of the switching probabilities and of the concurrence. 



4. Quantum noise 

To begin with, we consider the effect of quantum noise replacing in (|9]) Ea ~^ E^- 
The system dynamics is obtained by solving the Born-Markov master equation for the 
RDM. In the system eigenstate basis and performing the secular approximation (to be 
self-consistently checked) it takes the standard form [391 SO] : 

Pii{t) = -^Tim Pii{t) + ^ Tr,^i pmm{t) (14) 
kjit) = - {i^ij + f ij) p^j{t) . (15) 

The rates Fj^, Fjj and the frequency shifts Cjij — Uij, where depend 
respectively on the real and imaginary parts of the lesser and greater Green's functions 
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which describe emission (absorption) rates to (from) the reservoirs 

dte^-'{{E^{t)E^m) = -CEM--£EM (16) 

dte'-\{EMEc.m =i^CeA-^) + -SeA-^)- (17) 
In terms of the corresponding power spectra they read 



CeAuj) 



1 + exp {—(jj/kBT) 



oo 



^ / V „ / duj' Ce {^') , V 



OO 



where V denotes the principal value of the integral. Due to the symmetry of Ho + Hi, 
Eqs. ([3]) - (jlj), the only independent emission rates, Fj^, are Fio = r32, = Tsi, r2i, 
Fso, see figure [H Symmetric relations hold between the corresponding absorption rates 
Fjj. These processes originate from "effective transverse" noise, in particular transverse 
fluctuations (oc Xq.) enter the rates connecting the SWAP and the Z subspaces, whereas 
longitudinal fluctuations (oc Za) enter the intra-subspace rates, r2i, Fso, cfr eqs. (jlj) - (jTj). 
They read 

Tio = I (1 + sinv9) [Ca;^(a;io) + 6*^2(^10)] inter - subspace 

r2o = I (1 - sin v?) [Cx-i (w2o) + C'xa (^20)] inter - subspace , , 

Tso = \ sin^ (/? [C^j (W30) + (wgo)] intra - subspace 

r2i = \ [Czi(w2i) + C'zal^ai)] intra - subspace 

Absorption rates have the same form with CEa.^^im) replaced by CES~'^im)- The 
imaginary parts of the corresponding terms take similar forms. They enter the frequency 
shifts as reported in jAppendix C 



In the secular approximation, the SWAP and Z coherences decay exponentially 
with rates 

J- 12 = 7: [rio + Toi + + ro2 + ri2 + Fai] (21) 

f 30 = \ [rio + Toi + r2o + ro2 + Pso + Pos] + , (22) 

both inter-subspace and intra-subspace rates enter the decay of the SWAP and Z 
coherences. Note that the decay rate of the coherence puit) is only originated from 
dissipative processes (intra- or inter- effective transverse) since no pure dephasing 
processes (effective longitudinal noise) inside the SWAP subspace exist, cfr eq. (jlj). 
On the contrary, the coherences in the Z subspace also decay because of the effective 
longitudinal terms cos V5(|0)(0| -|3) (3 1 )(ii + 52), which originate P^* = I cos"^ ip[S, i{0) + 
5*22(0)]. This pure dephasing factor adds up to a decoherence rate due to intra-subspace 
effective transverse noise having the characteristic form [Tg^]"^ = (P30 + ro3)/2, as 
implied by (jHj), and to inter-doublet relaxation rates. 

Equations (jHj) for the populations do not decouple even in the secular limit. 
General solutions are quite cumbersome, so here we report expressions in the small 
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temperature limit with respect to the uncoupled qubits splittings, ksT <ti ^. In this 
regime, if the system is initially prepared in the state | H — ) = (|2) — |l))/-\/2, level 3 
is not populated, pssit) = 0, and the Z-coherences vanish, posit) = po3(0) = 0. The 
remaining populations are conveniently expressed in terms of the escape rates from levels 
1 and 2 

Tl = TlO + ; r2 ~ ^20 + ; (23) 

which enter the evolution of the populations in the following combinations 

r± = ^^^^ T \ v/(r? - Tlf + 4112121 . (24) 
For the chosen initial conditions the populations read 



2(r_ 

2(E 



k=± 



J- 12 

1 + 



6"^'=* (25) 



P22(t) = _^ . ^ [^12 - + n] e-"'=* (26) 
^' k=± 

Pooit) = 1 - (Pllit) + P22{t)) . (27) 



Note that, since ksT < <^ f^, thermal excitation processes internal to the SWAP 
subspace, expressed via the absorption rate ri2, cannot be neglected and r2i ~ ri2. On 
the contrary, inter-doublet thermal excitation processes are exponentially suppressed 
with respect to the corresponding decay rates, Fqi , ro2 <S Tio , r2o, with Fio ~ r2o- 
Thus the SWAP coherences decay rate (12T]) is approximately half the sum of the escape 
rates (El 



_ Pe _|_ re 

ri2 - . (28) 

In order to observe generation of entanglement inside the SWAP subspace it is necessary 
that relaxation processes to the ground state take place on a sufficiently long time 
scale. This is guaranteed when inter- and intra-subspace rates satisfy the condition 
riO)r2o ^ r2i5ri25 which requires that the spectra of the originally transverse and 
longitudinal fluctuations are Sx^{^) <^ Sz^{ujc) (from fl20l) ). In this regime, the SWAP 
decoherence rate ( l28i) is due to effective transverse processes internal to the subspace 

fi2~^[ri2 + r2i], (29) 

and the scales entering the populations take the approximate forms 

r+ ^ -[Fio + r2o] relaxation to the ground state (30) 

r_ ^ ri2 + r2i relaxation inside the swap subspace (31) 

with r+ ^ r_. Therefore, the time scales resulting from quantum noise, considering 
that ksT < < ^, when Sx^{^) < S'^^(wc), are 
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"Global" relaxation time to the ground state, analogous to the single qubit Ti. Its 
order of magnitude is the spectrum of transverse fluctuations in the computational 
basis at frequency Q 

8 4 
Tn = l/r+ ^ ^^^^^^ ^ ^^^^^^ = ^^^^^^ ^ ^^^^^^ , (32) 



where the approximate form comes from eq. fl20|) . where Uiq ~ a;2o ~ ^■ 

• Relaxation/decoherence times inside the SWAP subspace: They are due to 
"effective transverse" fluctuations inside this subspace, physically originated from 
longitudinal noise on each qubit at frequency Uc- Since there is no effective 
longitudinal noise in the SWAP subspace, relaxation and dephasing times are 
related by the typical relation T2^^^ ^ 2Tf^^^ where 

Tf'^^-P = l/r_ ^ l/(2fi2) (33) 

Ti^^' = i/r,2 - s.,{u,) + s.,{u,) ■ (''^ 

We now briefly comment on the validity of the secular approximation. It consists in 
separating the evolutions of elements Pijit) and pimif) provided that \uJij — uJim\ ^ 
where r denotes the typical evolution time scale of the system [39]. In the present case 
this condition is fulfilled if uj2i ^ T12. This constraint, on the other side, needs to 
be satisfied in order to observe generation of entanglement in the presence of quantum 
noise, as expressed for instance from anti-correlation of the probabilities (fTTl) and (fT2|) . 
Thus it can be regarded as a necessary condition, whose validity has to be checked case 
by case and requires (from fl2T]) ) 

S;,^ (Q) , (uc) < ujc . (35) 

In conclusion, entanglement generation in the presence of transverse and longitudinal 
quantum noise is guaranteed when 

5x.(f^) <5..(c^c) (36) 

T12 ^ ^ ' ^ ^ ' < uJc ■ (37) 

Under these conditions, the "global" relaxation time and the SWAP relaxation and 
decoherence times are given respectively by eqs. (132|1 and ( 133|) . ( 134|) . We note that, as a 
limiting case, the a/I — SWAP operation can be realized also releasing the condition 
r+ ^ r_, provided both rates are much smaller than the coupling strength Uc, 
and <^ ujc- For instance, when {ujc/^)Sx^{^) ^ Sz^{ujc) ^ uJc in Eq. 
iTio — r2o| + ri2 and the rates T± are still given by Eqs. ( l30l) . ( l3Ti) leading to 

Tr = l/r+ and Tf^^^ = l/r„. The SWAP dephasing time in this case also depends 
on trasverse fluctuations, T^^"^^ = l/fsi ^ l/[r+ + r_/2]. 
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Table 3. Noise characteristics deduced from single qubit data reported in ref. [3] . Since 
both qubits operate at qa.x = 1/2, they are presumably affected by similar polarization 
fluctuations, thus ^-^d Sx^ are the same as in table [2] Since the two qubits operate 
at different phase points, their phase spectral characteristics differ. From ref.|3] the 
power spectrum of phase fluctuations is Ss{uj) ~ [5.57 x fO^^J/w + 3.7 x fO^^'*[s] where 
the extrapolated crossover frequency is « 27rl0 MHz. Since 62 the process £2 oc AS2 
is Gaussian and characterized by 6*22(0;) = (iJj 2 sin ^2)^(^2,++ — ^2,--)^ Ss2{oj) 
(where $2,±± are defined in [Appendix A[ table IA2[) . Instead zi cx (AJi)^ is non 
Gaussian and its spectrum is ohmic, Szi{(j-!) = 2S^ / {ttH'^) uj coth{uj /{2KbT)), where 
5 = 1.6 X lO^s^^ In this table we fixed O w lO^^rad/s, ujc ~ 10~^^ and T = 40 mK. 





Qubit 1 


Qubit 2 


oV/ 


CT^, « 2 X iQ-'^n 


(7:^2 ~ 2 X iQ-'^n 




0-21 ~ 10-6 


0-22 « 6 X IQ-^fl 


si 


Sx^in) w4 X fO^s-i 


^^.(r?) w 4 X lo^s-i 


si 


5,,(c.,)«104 s-i 


5,2(wc)«5x 10^1 



/^.i. Switching probabilities 

In the secular approximation, for preparation at t = in | H — ) and for ksT <^ Q, the 
probabihties (fTTj) and (fT2il take the simpler form 

i'swi(t) = - ^ cos(^ [pn(t) + p22(t)] + Re[pi2(t)] + cos^ (^|) (38) 

P2{t) = - ^ cosv? [pn(t) + P22{t)] - Re[pi2(t)] + cos^ (^|) . (39) 
where pii(t) and p22{t) are given by eqs. ( l25l) . (l26l) and 

Pi2(t) = Pi2(0) exp {-t/Ti^^^^ - tuut} . 

Anti-correlation of the above probabilities directly follows from the coherence pi2{t) 
entering with different signs in Pswi(^) and in P2{t)- In order to check the efficiency of 
the gate we will therefore consider only the qubit 1 switching probability. Neglecting 



the frequency shift of LO21 (see [Appendix C[ ), -Pswi(^) can be approximated as 

Pswi(t) ^ I [e-*/^- - cos(a;2it) e"*/^^^"-^"] + cos^ (|) . (40) 

Here we explicitly see that, in order to perform the vi — SWAP operation, it is 
necessary that Tr , Tg^^^-f ^ 1/ujc- Efficient entanglement generation is guaranteed 
when Tr ^ j^swap ^ 1/uc, i.e. when transverse and longitudinal quantum noise 
are such that S^^^Q) <^ S^^^Uc)- Under this condition, the efficiency of the gate is 
limited by Tg^^^^, i. e. by longitudinal noise in the computational basis, Sz^{ujc), 
which is responsible for the short-times behavior. Decay towards the equilibrium value 
Pgwi(oo) = cos2((^/2) = (1 + cos(p)/2 = (1 - l/v^l + (cJc/2fi)2)/2 ^ {coc/^^f occurs 
in a time of the order of Tr, due to transverse noise in the computational basis. 
Capacitively coupled CPB-based qubits: The relevant rates (!20|) can be estimated in 



this specific case and are reported in Appendix A For the charge-phase two-ports 
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Figure 3. Left panel: Switching probability of qubit 1 for the noise levels and 
reported in table[3l The long-time decay is due to charge noise entering the populations 
of levels 1 and 2. Inset: the real part of the SWAP coherence which is responsible for 
the short-time behavior of the switching probability due to phase noise on qubit 2. 
Right panel: Switching probability of qubit 1 in the presence of charge noise with 
white spectrum, Sx-{uj) w 4 x lO^s"^. 



architecture, control is via the gate voltage, = CgVg/{2e) and magnetic flux dependent 
phase, 5, entering the Josepshon energy. The single qubit optimal point is at qx = 1/2, 
5 = 0. The resonant condition between the two qubits with a capacitive coupling is 
achieved by displacing one of the two qubits from the "double" optimal point. In order 
to limit the sensitivity to charge noise, resonance is achieved tuning 82 ~ 0.45. The 
noise characteristics are reported in table [31 where we note that S^S^) ^ Sz^{ujc). 
Because of the operating conditions, absorption (and emission) rates due to phase 
noise on qubit 2 dominate over rates due both to charge noise and to phase noise 
on qubit 1. Relaxation of the populations takes place on a scale 1/Tr = r+ ^ 
Fio ~ ~ 'S'xj(f2)/2 due to transverse (charge) noise on both qubits, whereas the 
SWAP coherence decay rate is dominated by longitudinal (phase) noise on qubit 2, 
IjTSWAP ^ f^^ ^ i^Y^2 + Tai) ^ T21 ^ S,,{io,)/4. The efliciency of the Vi - SWAP 
gate is mainly limited by phase noise on qubit 2, which gives T2^^^ ^ i/^2i ~ 100 ns. 
We note that, under these conditions, T.^'^'^^ is comparable with the decoherence time 
of qubit 2, T2 = l/[Sx^{n) /4: + S,2{0) /2] ^ 2/S,2{0) ^ 2/S,2{uJc)- The phase-dominated 
behavior is followed by a slower charge-dominated decay towards the equilibrium value 
-Pswi(oo) ~ {uc/4:QY. These features are illustrated in flgure [3] (left) where for 
comparison the real part of the SWAP coherence is shown in the inset. We note that 
the contribution of high-frequency charge noise to the efficiency of the y/i — SWAP 
gate is relatively small. Indeed if only polarization fluctuations were present, we could 
approximate 

Pswi(t) ^ ^ [1 - cos{co2,t)] e-*/^^^""" + cos^ (|) (41) 

and T2^^^ ^ 1/rio ~ 0.5 /xs. This situation is illustrated in flgure E] (right). 
Note that for the estimated noise flgures, the condition for the secular approximation, 
1^21 = ujc^ T12, is satisfled. 
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t(|is) 

Figure 4. Concurrence given by (|43)) for ujc/^ — 0.01 and for quantum noise values 
Sl and Sl reported in tabled Inset: At short times C{t) « 2|Im{pi2(i)}| (diamonds). 



4-2. Concurrence 



For the considered initial condition, the RDM takes the "X-form" 
concurrence is given by 



and the 



Cit) 



max 



0, V(pii(t) - P22(t))2 + (2Im[pi2(t)])2 - I sin(f\poo{t), 



sinv^lpooW - +P22(t))2 - (2Re[pi2(t)])2 



(42) 



At times shorter than the global relaxation time, Tr, the ground state is almost 
unpopulated, poo(^) ~ 0, and C(t) is given by the second term in (142|) 

1/2 



C{t) 



^e-2*/^« + sin2(a;i2t) e-^*/^^'""" 



sin<y9| (1 — e" 



(43) 



for t <^ Tfi, the concurrence is approximately given by the SWAP coherence C{t) ^ 
2|Im{pi2(t)}| = sin(a;i2t) e~*/"^2^'^'*^. Like the switching probabilities, the concurrence 
evolves with the SWAP coherence decay time, T2^^^, due to originally longitudinal 
noise. 

Capacitively coupled CPB-based qubits: The concurrence for the charge-phase 2-qubit 
gate is illustrated in figure HI We note that the long-time behavior is instead due to 
populations relaxation to the ground state and C{t) is given by the third term in (H2|) 

1/2 



C{t) ^ Isinc^l (1 



— I sm(p\poQ{t) — )■ I sin(y9| 



(44) 



The finite asymptotic value C{t — oo) ~ |siny9|, reflects the entangled thermalized 
state. Because of the interaction between the two qubits the phenomenon of 
entanglement sudden death does not take place 
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5. Adiabatic noise 

Let's consider now the effect of low frequency fluctuations, replacing in Ea — )■ 
E^ = Ea{t). In the adiabatic and longitudinal approximation [221 El] populations do 
not evolve and the system dynamics is related to instantaneous eigenvalues, Ui{E(t)), 
which depend on the noise realization, E{t). They enter the coherences in the eigenbasis 
of T-Lq in the form 

Pi.it) = P^j{0) [ V[Eis)] P[E{s)] e-^^ods^^Ms)) (45) 



where the probability of the realization E{s), P[E{s)], also depends on the measurement 
protocol. A standard approximation of the path-integral fHSl) consists in replacing 
Eaif) with statistically distributed values -E'a(O) = E^ at each repetition of the 
measurement protocol. The "static-path approximation" (SPA) [22] or " static- noise" 
approximation [13] gives the leading order effect of low-frequency fluctuations in 
repeated measurements. In the SPA, the level splittings Uij{E) are random variables. 



with standard deviation Sjj = (Sufj) — (Suij)'^, where 6uij = uJij{E) — Uij. The 
coherences (l45l) reduce to ordinary integrals 

p,,{t) ^ p,,(0) / dEP{E) e— -(^)* ^ p,,.(0) (e— -(^)*) , (46) 



where the probability density, in relevant cases, can be taken of Gaussian form, 
P{E) = U^P{E^) with P{E^) = exp[-El/2alJ/V2^aE^ [22]. The splittings uj,j{E) 
come both from "effective longitudinal" and from "effective transverse" terms in (jl]). 
This is analogous to a single qubit, where longitudinal noise gives the leading order 
linear terms and transverse noise is responsible for second order terms which dominate 
at the optimal point, where the first order longitudinal contributions vanish [T3| 1^ . 
Here, the Z-splitting ci;o3(-E') has a linear contribution due to the effective longitudinal 
noise (ii + Z2) cosv?(|0)(0| - |3)(3|) in dS}. The SWAP splitting U2i{E) instead, in 
the absence of leading-order effective-longitudinal intra-doublet terms in ([5]), comes 
from higher order contributions due to "effective transverse" noise. Evaluating them 
requires considering the complete Hilbert space of the coupled qubit system (inter- 
doublet processes included). The systematic approach to obtain these contributions 
consists in treating in perturbation theory effective transverse terms in Tii, where, in 
the adiabatic approximation, Xa and are replaced by classical stochastic fields Xa and 
Za- We obtain 

uj2i{xi,X2, zi, Z2) ^ujc- ^i^i + ^2) + ^(^1 ~ ^2)^ 

+ ^2)(^i + ^2) + - - ^2) (47) 
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/ !2 



UJo3{Xi,X2,Zi,Z2) ^ 2y ^ + fi2 _ cOS (^1 + Z2) 



+ ^ sin - cos [xl + ^2) + ^ sin <p {zi + Z2f , (48) 



where tamp = —uJc/{2Q). The SWAP-sphtting has been evaluated up to 4*'^ order |41j . 
whereas the Z sphtting expansion is considered up to second order since the fourth-order 
terms are much smaller, scaling with u^. The Z-splitting consists of a linear term due to 
the effective longitudinal noise, which dominates with respect to the quadratic term due 
to effective transverse intra-subspace noise (scaling with sirup). Transverse inter-doublet 
noise gives an additional second order contribution. 

Performing the average (H6l) (including in the S^'^ and 4*'^ order terms in ( 1471) only 
the contributions oc uj7^) we obtain the SWAP coherence in the SPA 



SPA 
Pl2 



where h{t) = {al^+al^ + iujc/t) {Vt^ / al + iUctY / {AVt^) and KqIH] is the K-Bessel function 
of order zero [l2]. In the absence of mixed terms in the expansion ( 148|) . different 
contributions to the Z-coherence factorize and lead to 

Pl^\t) = po3(0)eW^+-^* ; . . (50) 

l-^(Tfsin</^-cosy^) alt 

Here we assumed the same variance for the transverse noise components, = cFx2 = 
Instead we maintained and distinct, considering that the two qubits may operate 
at different working points (Appendix A). In (150|) the exponential factor comes from 
linear terms in fl48p due to effective longitudinal noise in the Z-subspace. The algebraic 
decay instead comes from quadratic terms in the expansion of cuos, due to effective 
transverse inter-doublet processes. 

We remark that the applicability of the Gaussian approximation to the fields depends 
on the relation between each and the fluctuations of the system's physical parameters. 
For instance, for a charge-phase yi — SWAP gate the resonant condition occurs for 
(5i = 0, ^2 7^ ( Appendix A ). Thus it is zi oc Elj{A6iY, where Si (not zi) is reasonably 
assumed Gaussian distributed. This would result in a modification of the terms oc 
in Eqs. (l49l) and (!50|) . For instance, the integral over zi in Pq^^ would reduce to 
{1 - iV2 cos^pazjt} "^/2, instead of exp { — | cos(y9|cr^^t2/2}. The quantitative effect in 
Pos^i^) (and similarly in pf2^{t)) is however negligible, because of the smallness of 
longitudinal noise at the optimal point, a^^ <^ a^^ (see table [3]). 

Validity regime of the SPA: The adiabatic approximation is tenable for times shorter 
than the relaxation times, for this problem this condition requires that t <^ j'SWAp^ 
The static approximation is exact for times smaller that I/7A/ (in case of a sharp high 
frequency cut-off). We verified that it is a good approximation also for times t > I/7M if 
7m < Q and the 1/f spectrum is originated from an ensemble of bistable fluctuators with 
switching rates G [7™,, 7m], leading to a 1//^ decay above 7a/ (numerical simulations 
and analytic first correction to the SPA [22]). 
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5.1. Minimization of defocusing: optimal coupling 

Exploiting the band structure of coupled nano-devices it is possible to reduce the 
influence of 1// fluctuations |13]. The basic idea of "optimal tuning" is to flx control 
parameters to values which minimize the variance of the splittings u!ij{E), Tj^-. This 
naturally results in a enhancement of the decay time of the corresponding coherence due 
to inhomogeneous broadening, i. e. of pff^it). This is simply understood considering 
the short times expansion of {exp {—i6uJij{E)t}) 

(e-'-'^)'>«l-,(^.,(i)>t-i(..,(E^^)t' (51) 
the short-times decay of the coherence in the SPA is therefore given by 



|(e-^'"'^(^)*)|^ Y^l-(S,,t)2, (52) 

resulting in reduced defocusing for minimal variance Sj^. For a single-qubit gate, 
the "optimal tuning" idea immediately leads to the well-known "magic point". In 
fact, if uJij{E) is monotonic in a region \Ea\ < 3cr£;^, we can approximate S^^- ^ 



2 

2 



^^l^^^o '^Ea^ thus the variance attains a minimum for vanishing differential 
dispersion. For the charge-phase two-port architecture, control is via gate voltage, q^, 
and magnetic- flux dependent phase S, thus Ea corresponds to the fluctuations Aq^, AEj 
and the optimal point, g^; = 1/2 , 5 = 0, is at the a saddle point of the energy bands [3]. 
When bands are non-monotonic in the control parameters, minimization of defocusing 
necessarily requires their tuning to values depending on the noise variances. For a multi- 
qubit gate, the optimal choice has to be done considering the most relevant coherence 
for the considered operation. Here we show how this program applies to the coherence 
in the SWAP subspace and partly to Po3^'^(^)- 

A key feature is that the SWAP splitting, 0721 in Eq. ( H7I) is non-monotonic in the 
small coupling Uc ^ ^. This is due to the fact that effective transverse fluctuations 
originate both from longitudinal and transverse noise, and give rise respectively to intra- 
SWAP and inter-doublet transitions, see ([5]), ([7]). For instance, second order corrections 
to ui from effective transverse intra-doublet processes are oc |(l|'Hi|2)p/(cJi — U2), 
from effective transverse inter-doublet transitions are ^j^^ 2 



Non-monotonicity in Uc results in a competition between 2"''^ and 4*^* order XQ,-terms 
in (H71) and in non-monotonic band structure, flgure [5] (panels (a) and (b)). Because 
of this subtle feature, identiflcation of the best operating condition necessarily requires 
consideration of the noise characteristics. Indeed the optimal coupling which minimizes 
the SWAP variance 

^4 



is given by 



- ^ { (I)' [i^. - ^l? + + aln^] + ^} (53) 



2 V a. 
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Figure 5. Dispersion of the SWAP and Z splittings for uic/^ = 0.02. Panels (a) and 
(b): 6tJ2i/^ from numerical diagonalization of + Hj. In panel (b) a zoom around 
the origin highlights the interplay of 2"'' and 4*'' order terms of the expansion (|47p . 
the barrier height is oc wj?. Panel (c): SWAP exact splitting (blue), expansion (|T7)) for 
X2 = 0, Zi = (dashed), 2"'' order expansion (dash-dotted green), Z splitting (red) 
from and single qubit dispersion (diamonds) . Panel (d) : Longitudinal dispersions 
in the SWAP (blue) and Z (red) subspaces. 



where we assumed ^ az^. The effectiveness of the optimal couphng choice has 
been discussed in details in ref. |13]. The advantage of operating at optimal coupling 
with respect to a generic Uc can be parametrized by the error of the gate at time 
te = 'it/2uc when system should be in the entangled state {ipe) = [| H — ) ~ i| — l")]/v^) 
e = 1 — {ipe\p(te)\4'e)- As showu in table m for two CPB-qubits, in the presence of 
moderate amplitude transverse (charge) noise, at the optimal coupling the error can be 
reduced even one order of magnitude with respect to a generic coupling. 

The splitting in the Z-subspace, on the contrary, is monotonic both in Xa and in z^, 
similarly to a single qubit operating at qa,x = 1/2 and at the two different 6a, figure 
panels (c) and (d). The most relevant contribution to Z- variance is given by the effective 
longitudinal noise 

^03 ~ cos^ ^ a^^ , (55) 
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Table 4. Error at the first \/i — SWAP operation at time for various couplings 
LUc/^ (numerical simulation of the coupled dynamics). The error Sa refers to a typical 
amplitude <7x — 0.02f2; the error Sb refers instead to a moderate amplitude = 0.0457. 
Longitudinal noise is here (Jz^ = 10~3J7. The optimal coupling. ([M)) . is in both cases 
at bJc ~ 0.0517. The error Sa is reduced by increasing ujc because of the comparatively 
large effect of phase noise. In the case of Sb, the most detrimental effect is from charge 
noise and the error is minimum at lUc- 
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Figure 6. Coherence in the SWAP-subspace (continuous red), in the Z-subspace (blue) 
and the single qubit-coherence (orange diamonds) for generic coupling lOc ~ 0.0117. 
The dashed red line is the coherence in the SWAP-subspace for optimal coupling, 
LUc ~ 0.0317. The Z coherence does not appreciably change by changing Wc and 
decays similarly to a single qubit at the double optimal point. Noise characteristics 
are reported in table |3l 

which can be reduced only by increasing the qubit 's couphng strength Uc, which is 
however hmited by single qubit splittings Q. Therefore, no special optimal point exists 
if the two-qubit operation involves the dynamics inside the Z-subspace. A comparison of 
the evolution of the coherences in the two subspaces and of the single qubit-coherence 
is shown in figure |6] for two CPB-based qubits. In order to keep the advantage of 
optimal tuning in the SWAP-subspace, the two-qubit operation should not involve the 
Z-subspace. 

Finally, we remark that the above results only follow considering the 4-level 
spectrum of the coupled devices. In fact, because of mixing between subspaces, limiting 
the analysis to noise projections inside the subspaces would miss a relevant part of the 
defocusing processes coming from originally transverse fluctuations oc Xa, see eqs. (jS]), 
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(jl]). This is clear if we consider the evolution of the SWAP coherences truncating 
the system Hilbert space to the SWAP subspace. Since the reduced Hamiltonian is 
([5]), effective transverse intra-doublet noise would lead to quadratic corrections to the 
SWAP-splitting [Jj]. In fact, treating in perturbation theory Za up to second order, we 
would get U2i{xi,X2, zi, Z2) ~ cUc + 2ir(^i ~ ^2)^ and the SWAP coherence would decay 
algebraically, as it is typically originated from " effective transverse" low frequency noise 

Pi2{t) OC / . (56) 

Under this approximation, reduction of defocusing could only be achieved by increasing 
the coupling strength. Similarly, reducing the analysis to the Z-subspace, cfr eq. 
would miss the quadratic dependence of the Z-splitting due to transverse noise. 

6. Time scales of the y/i — SWAP operation: interplay of quantum and 
adiabatic noise 

In the multi-stage elimination approach, where quantum noise is traced out first and 
adiabatic noise is retained as a classical stochastic drive, we have to replace in eqs. (l25ll 
- (1271) and in puit) = PuiO) Gxp {—t/T2^^^ — iuut}, Uij with ojij{E{t)) and perform 
the path-integral over the realizations of E{t). Note that the dependence on adiabatic 
noise enters also the rates, which should be averaged. If the dependence of the power 
spectra on the splittings is sufficiently smooth, it can be neglected and the effect of low 
frequency noise reduces to averaging phase factors in the coherences 0. Thus, provided 
that the dependence on E of and Fso can be neglected, the effects of quantum and 
adiabatic noise can be treated independently and lead to 

p,2(t)-pfr(t)e-/-^^'^^^ (57) 
Po3(t)^p^/'^(t)e-f3o*. (58) 

This result is confirmed by numerical solution of the stochastic Schrodinger equation 
for classical fluctuations leading to dynamical 1// noise (7^ = 27r x 1 Hz and 
7m = 2n X (10^ — 10^) Hz). The evolution of populations is only due to quantum 
noise (in the adiabatic approximation they do not evolve), thus they are given by Eqs. 
(j25l) - ( j27ll and the switching probabilities are given by Eqs. (j38l) . (j39l) where Pij{t) 
are replaced by (l57j) and (l58l) . These matrix elements enter also the concurrence which 
reads, for times t < T^, C{t) ^ V(pii(t) - P22(t))^ + (Im[2pi2(t)]) 2 - | sin (^| poo(t). 

The relevant time scales characterizing the efficiency of the Vi — SWAP operation 
depend both on the SWAP coherence and on the populations of the first three levels. 
Due to the interplay of adiabatic and quantum noise, the time dependence is not a 

II This is analogous to a single qubit in the presence of adiabatic transverse noise, which is equivalent 
to longitudinal quadratic noise j44] . 

% In the multistage approach, the classical variables Xa and Za enter parametrically also in the frequency 
shifts resulting from the solution of the Master Equation. Here we neglect this dependence, see 
[Appendix C[ 
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Table 5. Relaxation and decoherence times of the Vi — SWAP gate and responsible 
physical processes (long-to- intermediate time behavior of C{t) and of the switching 
probabilities) . 



Time scclIg 


Physical origin 


RelcLXcitiorL to the ground stcitc 


Transverse noise 


Tr « i/r+ 


at frequency f2 


SWAP relaxation time 


Longitudinal noise 


rpSWAP ~ x/p 


at frequency U!2i 


SWAP decoherence time 


Longitudinal noise 




at frequency uj2i 


SWAP dephasing time T^^^^* 


Low frequency noise 




(longitudinal AND transverse) 


SWAP total decoherence time 


if Tn » i. e. 


rpSWAP _ ^i/2Tf^"^^ + i/xSWAP*-^-i 


s^in) « 5,(w2i) 



superposition of exponentials. We can distinguish two time regions: a asymptotic long- 
time regime and a intermediate-to-short time regime. 

The asymptotic behavior is entirely due to populations relaxation to the ground 
state, it is exponential and takes place with the "global" relaxation time Tr ^ l/^-i-, 
resulting from "effective transverse" inter-doublet processes, whose order of magnitude 
is the spectrum of transverse fluctuations at frequencies of order Q. In order to avoid 
leakage from the SWAP-subspace, any two-qubit operation has to take place on a time 
much shorter than Tr. This constraint also applies to the SWAP decoherence times. 
The intermediate-to-short time behavior, t <^ Tr, gives more relevant information 
on the gate performance. In this time regime, relevant quantities are populations 
and coherences in the SWAP subspace. We distinguish a intermediate time regime, 
characterized by Tf^^P ^ l/r_ and T^^^p = 2Tf^^^ ^ 1/fu due to effective 
transverse quantum noise inside this subspace, physically originated from longitudinal 
noise on each qubit. Adiabatic noise leads to additional decay of the SWAP coherence, 
Pi2{t) ~ Pi2^{i) e'^l^'^'^^^ . The resulting defocusing is analogous to a "pure dephasing" 
process, we may name the typical time scale j'^^^^^*^ defined as the time at which 
|p^f^(T2^^^^*)| = e-^ The time T^'^^p* is found by numerical inversion of fli9|) . 
The intermediate-time behavior of the two-qubit gate is characterized by \jT^'^^P and 
XjT^'^AP = l/(2Tf^^^) + l/T^^^P*, analogously to a two-state system. These time 
scales and the responsible processes are summarized in table [5l 

On the other side, in a quantum information perspective, it is important to estimate 
the time behavior (of concurrence or of switching probabilities) at times of the order of 
te = tt/2uJci the first moment in which the entangled state is reached by free evolution. 
In order to generate entanglement within the SWAP subspace it is necessary that 
tf. <^ T2^^^. It is therefore relevant to estimate the system behavior at short-times 
t <^ T2^^^. Both for the concurrence and the switching probability, the SWAP- 
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t (|a,s) t (|j.s) 

Figure 7. Qubit 1 switching probability Pswi(^) (red) and probability P2{t) (blue) 
to find qubit 2 in the initial state |— ) in the presence of 1// and white noise for 
Wc/il = 0.01 and for optimal coupling Cuc = O.OSfi (gray). To emphasize the robustness 
of the optimal coupling choice here we consider 1// charge noise of large amplitude, 
(Tj; = 8 X lO^^O. The remaining noise figures are those reported in tabled Left panel: 
Plot of Pswi(i) showing the exponential short-time behavior at a)c and the algebraic 
decay for generic coupling. Right panel: Pswi{t) and P2 anti-phase oscillations for Hjc 
(main), Uc/i^ = 0.01 (inset). 



coherence rules the short-time hmit which is approximately 

IP12I oc 1 - ^y/^^p - ^ ^21^^ • (59) 

The leading contribution, linear or quadratic, does not only depend on the noise 
amplitude but also on the operating point. In fact Uc enters the SWAP-splitting variance 
E21 and, in principle, also the SWAP decoherence time due to quantum noise 2Tf^^^. 
This is expected to be a smooth dependence. For the present considerations we assume 
a white spectrum in the relevant frequency range. For optimal coupling the short-time 
behavior is linear (since the effect of low frequency noise is considerably reduced) whereas 
for a different coupling strength the behavior is typically algebraic (see figure [7]). 

The short-time expansion is valid at tg if te <C T2^^^. In figure M we consider 
two illustrative cases. For the expected noise characteristics as reported in table [3l for 
couplings n/100 <ujc< ^/lO we have Ssite < te/2Tf^^^ < 1. Thus the short-time 
expansion is valid up to tg. Note that for the optimal coupling l/2Tf^^^ is about one 
order of magnitude larger than S21 (figure |8] left panel), thus IP12I ocl-t/2Tf^^^. For 
a larger amplitude of 1// transverse noise, ax ~ O.OSfi, we have te/2Tf'^^^ < S2ite ^ 1 
(figure [8] right panel). Also in this case the short-time expansion holds up to te- In this 
case however even for optimal coupling the linear and quadratic terms are comparable. 
The various possible short-times behaviors and validity regimes are summarized in 
table [6l 
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Figure 8. Plot of Ei2(wc)ie (red) and of ie/(2rf ^^^) (blue) as a function of Wc/^^ 
(logarithmic scale) for cr^/ri = 0.02 (left panel) and ax/^ — 0.08 (right panel). The 
dashed line marks the optimal coupling Eq. ([54| . 



Table 6. Short times behavior t ^ j^swap C{t) and of the switching probabilities 
and validity conditions of the expansion up to tg. 

Short times Validity Expansion 

expansion conditions up to if 

|p2i(0l ~ 1 - 5rprxp Optimal coupling and S21 • 2rf^^^ < 1 ^j^^tp < 1 
\P2i{t)\ ~ 1 - Generic coupling or E21 • 2Tf '^'^^ > 1 < 1 



7. Conclusions 

In the present article we have identified the relevant decoherence times of an entangling 
operation due to the different decoherence channels originated from broad-band and non- 
monotonic noise affecting independently two qubits. Results depend on the interplay of 
noise at low and at high frequencies (with respect to both the single qubit splittings, Q, 
and the qubits coupling strength, Uc) and on the symmetries of the qubit-environment 
coupling Hamiltonian. In addition, "optimal" operating conditions for the universal 
two-qubit gate have been identified. Even if the relevant dynamics is within the SWAP 
subspace, the important time scales (both due to quantum and to adiabatic noise) can 
only be predicted considering the whole multilevel nature of the coupled systems. 

In particular, relaxation processes from the SWAP subspace to the ground state only 
depend on transverse noise at frequency Q. Decoherence processes internal to the SWAP 
subspace are instead originated from longitudinal noise at frequency Uc- This apparently 
counter-intuitive result simply follows from the symmetries of the Hamiltonian expressed 
in the eigenbasis of coupled qubits, cfr eqs. (E]) - (I7j). As a consequence, the conditions to 
generate entangled states within the SWAP-subspace involve the spectra of transverse 
and longitudinal fiuctuations at two different frequencies. An efficient ^/i — SWAP 
operation can be realized when S^ai^) ^ Szai^c)- The long-to- intermediate time 
behavior and the relevant decoherence times are summarized in table |5l 
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Defocusing processes instead originate both from transverse and longitudinal 
adiabatic noise. Remarkably, consideration of the coupled systems band structure allows 
identification of operating conditions of reduced sensitivity to low frequency fluctuations. 
An "optimal coupling" can be identified where defocusing is minimized. The error at 
the first \/i — SWAP operation at optimal coupling can be reduced of a factor of 4-to-lO 
with respect to operating at generic coupling strength. The possibility of optimal tuning 
is ultimately due to the absence of effective longitudinal noise in the SWAP-subspace 
and to the interplay of effective transverse intra- and inter-doublet fluctuations. The 
orthogonal Z-subspace instead, because of the presence of effective longitudinal noise, 
turns out to be much more sensitive to low-frequency noise which cannot be limited by 
properly choosing system parameters. Therefore, populating this subspace may severely 
limit the gate efficiency. 

The short-times behavior of the relevant dynamical quantities crucially depends 
on adiabatic noise and its interplay with quantum noise. The dependence turns from 
quadratic to linear depending on the largest component, parametrized by the SWAP- 
splitting variance resulting from adiabatic noise, S21, and the SWAP-decoherence rate 
due to quantum noise [2Tf'^^^]~^, as summarized in table El 

The considerable protection from 1/f noise achievable in the SWAP-subspace for 
optimal coupling may suggest to use this subspace to encode a single quantum bit, 
similarly to a decoherence- free-subspace |15]. Such an encoding would be convenient 
if the SWAP decoherence time was about two orders of magnitude larger than the 
qubit decoherence time T2. In fact, because of the smaller SWAP oscillation frequency, 
Uc ~ 10~^*-"^''f2, the quality factor of a single qubit rotation within this subspace would 
be Q = T2^^^V ^ 10-i(-2)T2^^^^(], therefore Q > T2Q if Ts^^^^ > 10'^^^T2. 
Realizing this condition requires improving the relaxation times with respect to present 
day experiments. A more realistic possibility is instead to employ the SWAP subspace 
as a single qubit quantum memory. The requirement in this case would be satisfying 
the weaker condition, Tg^^^-P > 7^2. 

Finally, we would like to comment on the effects of selected impurities strongly 
coupled to the two-qubit gate. The detrimental effect of charged bistable fluctuators 
on single qubit charge or charge-phase gates has been observed in experiments [21 |3] 
and explained in theory [221 iSE]. Effects on a two-qubit gate may even be worse, as 
reported in the recent experiment on two coupled quantronium [38]. Explanation of 
the rich physics which comes out when selected impurities couple to the nano-device 
is beyond the scope of the present article. In the following Appendix B we illustrate 
the possible scenario when an impurity considerably changes the 4-level spectrum of the 
coupled qubits and we speculate on the possibility to limit its effects by proper tuning 
the system parameters, somehow extending the optimal tuning recipe. The usefulness 
of such a choice however critically depends on the interplay with quantum noise at 
intermediate frequencies, an information unavailable in present day experiments. The 
analysis in Appendix B is the first step towards the identification of parameter regimes 
where a multilevel nano-device may be protected also from strongly coupled degrees of 
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freedom of a structured bath. 
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Appendix A. Capacitively coupled Cooper-Pair-Box based qubits 

The Hamiltonian T-Lq in Eq. ([1]) is the central model of any non trivial two-qubit gate 
based on a fixed coupling scheme. In particular, it describes coupled superconducting 
qubits in the various implementations [5]. With Tii, given in ([2]), it includes independent 
fluctuations responsible for incoherent processes of different physical origin. In this 
Section we derive T-Lq + "Hi for capacitive coupled Cooper-Pair-Box-based (CPB) nano- 
devices. The CPB is the main building block of many superconducting qubits |46] . 
Here we consider two charge-phase qubits (quantronia) [3] electrostatically coupled via 
a fixed capacitor, as schematically illustrated in figure lAll 

Each qubit is operated via two control parameters, the gate voltage Vg and the 
magnetic flux across the junction loop, They enter the CPB Hamiltonian via the 
dimensionless parameters = CgVg/{2e) and S = n^x/^o- 

-HcPB = E'ciq- q.f-Ej{5) cos ^, (A.l) 

here = 2e^/CE and Ej{6) = EjCos6 is the Josephson energy, modulated by the 
phase 6 around the zero phase value Ej. Charge q and phase (p are conjugate operators, 
[(f, q] = i. The loop capacitance is the sum of the gate Cg and the junctions Cj 
capacitances. The two CPBs are coupled by inserting a fixed capacitance Cq between 
the two islands. The presence of the coupling capacitance leads to a renormalization of 
the CPB charging energies E^^c = E^q{1 — Ct/Cq,_s), being E^^ the CPB charging 
energy of the uncoupled qubit a, 1/Ct = I/Cq + l/C'i,^ + 1/C*2,s the total inverse 
capacitance of the device. The coupling energy is Eqc = (2e)^CT/(Ci^sC'2,s), and the 
full device Hamiltonian reads 

no = HcFBl ® I2 + II ® HcPB2 + He (A.2) 
"HcPBq = Ea,c[qa " qa,^laf " -EQ,j(5a) COS<^«, 

Uc = Ecciqi - qi,M ® - g2,xl2], (A.3) 

where the control parameters qa,x and 6a, and the charge qa and phase operators (pa 
play the same role as in the single qubit case. At sufficiently low temperatures, the 
evolution of each CPB approximately takes place within the bi-dimensional subspace 
spanned by the lowest energy eigenstates of "HcPBa, with a splitting depending 
on the control parameters [17], "Hq = Va'HcPBa'Pa = —^^aO'az, where the projection 
operator reads Va = a\+){+\a + Q|-)(-|a and aaz = a|-)(-|Q - a\+){+\a- The 
restricted dynamics of the coupled CPBs can be described in a pseudo-spin formalism by 
projection in the eigenstates basis u) = ® |z^2)}- In this subspace the charge and 
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Readout 1 




Readout 2 



Figure Al. Electrical circuit of two quantronia, labeled as 1 and 2, coupled 
via a fixed capacitor Cc- Each qubit has the characteristic two-port design, with 
control parameters the gate voltage Vg^a and the magnetic flux $2^_q, threading the 
superconducting loop [Sj. 



Josephson operators are expressed in terms both of a^z and of the transverse component 

(yax= Qj+)(-U + a|-)(+U as follows 



P„(C0S 0a)'Pa 

where 



1 



-a, 



ax I 2 ("^"i 



The restriction of the device Hamiltonian to z/)} is El 
(Pi ® V2)Hy,{Vi ® Pa) = ^ ^ 



+ (A.5) 



E, 



cc 



E 

^CC JJ^ 



;{qa,++ - Qa, )o-az + ^q.-H-O-q 



(A.6) 



Note that because of the coupling, in addition to interaction terms between the two 
qubits, each qubit is effectively displaced from its own operating point (second term in 
(lAj])). 

Because of fluctuations of control parameters the device couples to the external 
environment. The form of interaction terms can be easily deduced considering classical 
fluctuations of qa,x and Sa, i.e. by replacing in f lA.2p and flA.Sp . qa,x — > Qa,x + ^?q,x and 
^ + ^^a- Wc therefore obtain "Hd — + SUd, where 

5Hd = SncPBl ® I2 + II ® 5HcPB2 + SKc- 
Wc have used the shorthand notation ctq to indicate ctq, ® I^, a ^ (3. 
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Fluctuations of each CPB Hamiltonian and of the couphng term read 

= XaQa + COS 0a 

S'Hc = 92X2 gi ® II2 + fi-l^lll ® q2 , 



(A.7) 
(A.8) 



where is related to gate charge fluctuations, Xa = —2Ea^c^<la,^ and to 
fluctuations of the phase or equivalently of the Josephson energy, Za = —AEa^j, where 
AEa^j = Ea,j{Sa + A6a) — Ea^j{Sa) , and we put 9a = EccI {'^Eafi). In the computational 
subspace the additional terms reduce to 

5T-ia = "2"^"''^"'++ ~ )^az + Xaqa,+-0'ax + 2^"(^"'++ ~^ 

-^Zai<^a,++ - '^a,^^Wz + Za<^a,+^(rax + ^^a($a,++ + $a, — )Ia (A.9) 



and 



SHc = 92X2 



+9iXih 



-^(92,++ - 



91, 



® I2 



Q2, 



)0-2z + g2,+-0-2a 



2(92,++ 



+ g2,— )l2 



. (A.io) 



In conclusion, the coupled CPBs Hamiltonian at a general working point, in the four 
dimensional subspace z/)}, including fluctuations of the control parameters and 
neglecting constant terms, takes the form "Hd + SHb with "Hd given in eq. ( 1A.6I) and 

'1 



a, 



a,++ 



))aaz + {Xaqa,+- + ^«$a, 



-92X2 



jO-lz + gi,+-Cri:i 



+ ^71^1 



g2,++jo"2^ + q2, 



-0-22 



(A.ll) 



Note that due to the capacitive coupling a cross-talk effect takes place, gate charge 
fluctuations of qubit a being responsible for fluctuations of the polarization of qubit 
(3. Effects are scaled with the coupling energy Ecc, thus they are expected to be less 
relevant compared to fluctuations acting directly on each qubit. A detailed analysis of 
cross-talk and correlations between gate charge fluctuations has been reported in [33] . 
In our analysis we disregarded cross-talk and correlations between noise sources acting 
on each qubit. 



Appendix A.l. Choice of the working point 

The implementation of a two qubit gate in a fixed coupling scheme requires the qubits 
to be at resonance during gate operation. Because of experimental tolerances on bare 
parameters (about 10%), the resonant condition can only be achieved by a proper choice 
of the single-qubit working points via tuning ^ and 5a- This is a critical choice since 
at least one qubit has to be moved away from the working point of minimal sensitivity 
to parameters variations, the "optimal point", = 1/2, (5a = [3]. Since CPB-based 
devices are severely affected by low-frequency charge noise it is desirable to maintain 
the charge optimal point, g^^x = 1/2 and modulate resonance/detuning by tuning the 
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Table Al. Typical values of the CPBs parameters [H^. is the energy splitting of 
the single qubit. 



Parameter 


Qubit 1 


Qubit 2 


£;0 (GHz) 


9.997 


10.41 


E'j (GHz) 


14.16 


15.62 


Cs (fF) 


7.73 


7.42 


n° (GHz) 


12.666 


13.81 


Ea^c (GHz) 


9.92 


10.33 


n (GHz) 


12.645 


13.79 



Table A2. Matrix elements of the charge and Josephson operators in the two lowest 
eigenstates basis of each CPB Hamiltonian. The physical parameters of each CPB are 
given m table EH 

Qubit 1 Qubit 2 

qa, + + - Qa, 

qa,+- -0.5962161 -0.5891346 

!>„,++ - 0.7049074 0.7396952 

+ - 



phases 6a- The most reasonable choice consists in maintaining one qubit at its own 
phase optimal point, 6i = and tune the phase of the other qubit, 62. For set of 
parameters close to those planned in experiments [3H] (see table lAll) resonance occurs 

for 62 ~ 0.455. Since at g^^x = 1/2 it results qa,++ = Qa, = 1/2 and = 0, the 

projected charge ( 1A.4I) is transverse and the Josephson operator ( ]A.5|) is longitudinal 
with respect to each Hamiltonian "Hq,. Therefore, the truncated Hamiltonian is of the 
form "Ho + given by eqs. ([1]) and ([2]) respectively 

"^0 = ® I2 - II ® + Ecc qi,+~ q2,+- Oxx ® (^2x (A.12) 



5no = Yl 



-Za{^a, ^a,++)(^az + Xaqa,+-(yo 



(A.13) 



where qa,+-i etc. denote the specific values taken by these matrix elements 

at the resonance point obtained for a specific choice of 82- Values of the charge 
and phase matrix elements at this working point are reported in table IA2I and lead 
to Eqc = 0.18 GHz. The interaction Hamiltonian Hi results from STio considering 
the quantized version of the classical fluctuations Xa = —4:Ea,cQa,+-^(la,x and Za = 

{^a,++ - ^a, )^Ea,j. Where, since 5i = we have AEij = -0.5£^5^ j(A5i)^ whereas 

for 7^ it results AE2,j = -E° j(sin(52) A62. 

The accuracy of the truncation of the multistate device Hamiltonian to the lowest 
four eigenstates has been checked numerically by exact diagonalization of ( ]A.2|) using 
the numerical values of the physical parameters reported in table lAll and considering 
14 charge states for each qubit, see figure EH 
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Figure A2. The five lowest eigenenergies of the coupled two CPB Hamiltonian (jA.2[) 
at charge optimal points <ji_x ~ 'Z2.x = 1/2 versus the qubit 2 phase 62, being Si = 0. 
The dashed black and red lines refer to the uncoupled and coupled case, respectively. 
The 5-th eigenvalue lies sufficiently far in energy to be neglected in our analysis. In the 
inset the degeneracy point is shown with the small anti-crossing due to the coupling 
term at 62 = 0.455. 




Appendix A. 2. Structured noise in CPB-based circuits 

Fluctuations of the control parameters go^x, have different physical origin. They partly 
stem from microscopic noise sources, partly from the circuitry itself. The resulting 
stochastic processes are sometimes non-Gaussian, therefore complete characterization 
of the processes requires in principle higher-order cumulants. Such a complete 
description is often unavailable in experiments, whereas it is usually possible to provide 
a characterization of the power spectrum of the various stochastic processes. Noise 
characteristics for the two quantronia are summarized in table [3] (section Hj). 

Experiments on various Josephson implementations have revealed the presence 
of spurious resonances in the spectrum which manifest themselves as beatings in 
time resolved measurements [3]. In charge and charge-phase qubits they are due to 
strongly coupled (SC) background charges and may severely limit the reliability of these 



devices [351 ESj . In the following Appendix B we will discuss the detrimental effect of a 



single strongly-coupled impurity on the visibility of a \Ji — SWAP operation. 

Finally, we mention that in the present analysis we have not explicitly considered 
fluctuations of the Josephson energy due to critical current fluctuations. The effect of 
noise may be straightforwardly included in our analysis within the multistage approach. 

Appendix B. Effect of a strongly coupled impurity 

Experiments on various Josephson implementations have revealed the presence of 
spurious resonances in the spectrum which manifest themselves as beatings in time 
resolved measurements [H [15], [161 [21] • In charge and charge-phase qubits they are due 
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to strongly coupled (SC) background charges and may severely limit the reliability of 
these devices [35], [36]. Several experiments have shown that these impurities can be 
modeled as bistable fluctuators (BF), switching between states and 1 with a rate 
7. Often the qubit-BF coupling is transverse, i.e. of the form —{v/2)axC,(t) with 
^(t) = {0,1} [IE]. Each bistable state corresponds to a different qubit splitting, fl 
for = and fl' ~ Q[l + 0.5(i;/fi)^] for ,^ = 1 (t; <^ Q). The parameter quantifying 
the strength of the qubit-BF coupling is gi = {Q' — Q)/'~f [35]. When gi ^ 1 the 
impurity is strongly coupled and it is visible both in spectroscopy and in the time- 
resolved dynamics. Similarly, a single BF acting on one of the two qubits forming a 
universal gate induces a bi-stability in the SWAP-splitting: 0021 = Uc foi =0 and 
u'21 ~ Wc[l — 0.5(t'/f2)^ + v'^ / [Su'^Q'^)] for ,^ = 1 (in this case the stronger condition on 
the qubit-BF coupling has been assumed v <^ Uc). The strong coupling condition for the 
a/I — SWAP gate is therefore g2 = (a;2i — to'2i)/7 ^ 1- Considering the above expansions 
and the constraint Q > Uc the two ratios satisfy the condition gi > g2- Therefore the 
visibility of a specific BF is expected to be reduced in a two-qubit operation with respect 
to a single qubit gate. 

The effect of a single BF fluctuator depends on the way experimental data are 
collected. For spectroscopic measurements the single-shot time is about 200 ns and the 
acquisition time of a single point in the spectrum requires about 4 x 10^ repetitions, 
resulting in a recording time for a single data point of about 10~^ s [3H]. Therefore 
a single BF will be visible in spectroscopy if its switching time is 0.5 x 10~'^s < 
1/7 < 2 X 10~^s. For instance a BF switching at ~ 10 kHz can be averaged during 
the spectroscopic measurements and it is visible if it is sufficiently strongly coupled 
to one qubit, for instance this is the case for v/Q = 0.1 > ujc/^- The expected 
effect of this impurity in spectroscopy is shown in figure IBll where it results in the 
simultaneous presence of two avoided crossings depending on the impurity state. The 
occurrence of a similar BF represents a major problem for charge and charge-phase 
implementations [38]. Effects are also visible in time resolved measurements. If the 
considered BF is one out of the several ones responsible for 1// noise, the global effect 
of the structured bath is a considerable reduction of the oscillation amplitude, as shown 
in figure IBll Under these conditions, in principle, an "optimal coupling" can still be 
identified. Since the main problem is the beating pattern, optimal tuning is defined 
by the condition that the average {uj'21 — UJ21) vanishes, rather than minimizing the 
SWAP-splitting variance. In the presence of only charge noise, this condition leads to a 
modified "optimal coupling", Cjc = y^ct^ + at which effectively the SWAP visibility 
is improved (figure iBTl) . The efficiency of such a choice however critically depends on 
the details (presently not available) of transverse noise at frequencies in the kHz - MHz 
range. Relaxation processes due to quantum noise at these frequencies may in fact 
represent an additional liming factor for the gate fidelity. 
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Figure Bl. Left: Simulated spectroscopy of the SWAP splitting UJ21 for qubit-qubit 
coupling uJc/^ = 0.01 (blue) and for uncoupled qubits (black). In the presence of 
a BF coupled to one qubit with strength v/fl = 0.1 the coupled qubits {oJc/^ = 
0.01) eigenenergies are shifted (red). Right: Numerical simulations of the switching 
probabilities in the presence of the above BF (orange/cyan), and of the BF plus 1// 
noise with ax = 0.02r2 (blue and red). For optimal coupling, ojc = -v^ct^ + v'^/A, a 
considerable recover of the signal can be obtained (gray). 



Appendix C. Frequency shifts 



In this Appendix we report the frequency shifts entering the off-diagonal elements of the 
RDM as resulting from the solution of the Master Equation in the secular approximation 
reported in section |H The frequency shifts take the simple form 



^03 



W03 



(C.l) 



For the SWAP coherence the different contributions read 
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£12 = 




(^12)] 



From single qubit measurements, as reported in section Appendix A. 2 , charge 
noise at frequencies of order of Vt is white. Assuming that also the phase 
variables A(5a have white spectrum at frequencies of order cjc, results in Sz^ioj) = 



{E%sm62rSs,{uj) 
Szi (uj) 



2q^ f, ,\ ^ g ^ ]^g7g i_ rjij^g spectrum of Zi takes instead the ohmic form 
2SV(7rfi2)a; coth(u;/(2i^BT)), where ^ = 1.6 x lOV^ Under these 
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conditions, phase shifts due to polarization noise and phase fluctuations on qubit 2 
identically vanish (from f|T9|) ). A frequency shift contribution results from ohmic phase 
noise on qubit 1, similarly to a qubit affected by transverse noise [23], as it can be argued 
from ([5]) 

8 4:7!- 1 + e '^/^Bi _ 

The shift depends on the high-frequency cut-off of the spectrum and on the temperature. 
Explicit forms are known in the literature, see [23]. The damping strength parameter 
entering the ohmic spectrum, commonly denoted as K, is in the present case K = 
5'V(7rfi^) ^ 10"^. The SWAP-splitting shift scales with K and it will therefore be 
extremely small. Possible logarithmic divergences with the high-frequency cut-off are 
not included in this analysis since there's presently no indication on the features of the 
spectrum at those frequencies. For this reason frequency shifts have been neglected both 
in Section H] and in Section O 
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